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formula) , gauge transformations of the wave function and generalized symmetries of the 
associated integrable ODE. The paper contains a detailed explanation of the immersion 
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1. Introduction 

The language of group theoretical methods for investigating soliton surfaces associated 
with integrable models is a very useful and adequate tool for studying the main 
features of various problems appearing in diverse areas of theoretical physics. The 
possibility of finding an increasing number of 2D surfaces immersed in Lie algebras 
can be discovered after analysis of the group properties of the models. The main 
advantage of such an approach appears when group analysis makes it possible to 
construct algorithms for finding particular classes of 2D surfaces without referring to any 
additional considerations but proceeding instead directly from the given model under 
consideration. A systematic method for constructing surfaces and their continuous 
deformations under various types of dynamics has been extensively developed by many 
authors (see e.g. [91 QUI [12], [131 E])- A broad review of recent developments in this 
subject can be found in several books [H |2j HJ [7J [TTj, (20J |2l] and references therein. 

The methodological approach assumed in this work is based on the Fokas-Gel'fand 
approach, as presented in [14] . providing a symmetry characterization of deformations 
of soliton surfaces using the formalism of generalized vector fields and their prolongation 
structure. The identification of necessary and sufficient conditions for the existence of 
2D surfaces in terms of invariance conditions for generalized symmetries allows us to 
integrate the immersion functions explicitly in terms of conformal transformations in the 
spectral parameter of the linear spectral problem (LSP), gauge transformation of the 
wave functions and generalized symmetries of the model and its linear spectral problem. 
The results obtained for systems of PDE's were so promising that it is worthwhile to 
try to adapt this method and check it effectiveness for the case of ordinary differential 
equations (ODE's) associated with elliptic functions. This is, in short, the aim of the 
present paper. 

The crux of the matter is that we consider an ODE in the dependent function u 
and the independent variable x which can be written in the matrix Lax representation 
[T5| 123 ES] involving the differentiation of the Lax pair with respect to x only. This type 
of Lax pair of equations is obtained by a compatibility condition of some LSP for which 
an auxiliary variable, say y, has been introduced in the wave function $. In this setting, 
the wave function becomes a function of A, x, y, u and the derivatives of u with respect 
to x. The main advantage of this procedure is that, in using the invariance criterion for 
the generalized symmetries in terms of their prolongation structure, it leads to simple 
formulae and allows us to write the explicit form of the soliton surfaces. 

The plan of the paper is as follows. Section [2] contains a brief account of 
basic definitions and properties concerning the Fokas-Gel'fand formula for immersion 
functions in Lie algebras. We concentrate on the study of soliton surfaces associated 
with ODE's in Lax representation for elliptic functions. In section [3j we investigate 
in detail the Lax pair and find the general solutions for its wave functions for these 
types of second-order ODE's. These results are then used in section [4] to formulate 
soliton surfaces associated with conformal symmetry in the spectral parameter, gauge 
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transformations of the wave function, and generalized symmetries of the associated 
integrable ODE and its LSP. In section [5] and section [6j we present examples of 
applications of our approach to the case of the Jacobian and P-Weierstrass elliptic 
functions respectively and calculate their geometric characteristics. Section [7] contains 
final remarks concerning soliton surfaces associated with the elliptic functions, identifies 
some open questions on the subject and proposes some future developments. 

2. Application of the Fokas-Gel'fand formula to ODE's in Lax 
representation 

In this paper, we construct soliton surfaces immersed in Lie algebras using the Fokas- 
Gel'fand formula for immersion [13], as formulated in [TJ], applied to ODE's. For this 
purpose, consider an ODE 

A[u] = A(x,u x ,u xx ,. . .) = (1) 

which admits a Lax pair with potential matrices L(X, [u]), M(A, [u]) taking values in a 
Lie algebra g which satisfy 

D X M + [M, L] = 0, whenever A[u] = 0. (2) 

Here, A is the spectral parameter. In what follows, we make use of the prolongation 
structure of vector fields as presented in the book by P. J. Olver [19j. For derivatives of 
u we use the standard notation 

D x u = u x , D x uj = uj )X , J = (x, x) 

and, for functions depending on the independent variable x, dependent variable u and 
its derivatives, the following notation has been used 

f[u] = f(x, Ux, u xx , ...), /(A, [it]) = /(A, x, u x , u xx , ...). 

The total derivative in the direction of x takes the form 

d d 

D x = D x + u x — + u xx - h (3) 

ou au x 

This Lax pair equation (|2| can be regarded as the compatibility conditions of a 
linear spectral problem (LSP) for wave functions $ taking values in the Lie group G 
with independent variables x and an auxiliary variable y, and spectral parameter A. For 
the purpose of symmetry analysis, we allow <3> to also depend on the dependent variable 
u as well as its derivatives with respect to x. The LSP is written 

D x $(\,y, [u]) = L(X, [u))$(\,y, [«]), D y $(\,y, [u]) = M(A, [«])$(A,y, [«]). (4) 

Here, since the function u is independent of y, the total derivative in the direction y is 
given simply by 

°< = k 

Note that, because L and M also do not depend on the auxiliary variable y, the 
compatibility conditions for ^ are of Lax form ([2]). 
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With such an LSP, we can then apply the results of [14] . to construct soliton surfaces 
immersed in the Lie algebra g. That is, suppose that L, M G q and $ G G satisfy the 
LSP Q and its compatibility condition ([2]), then there exists a g-valued function F 
with tangent vectors given by [13J 

D X F = D y F = Q^B® (5) 

for any g-valued functions A(\,y, [u]) and B(X,y, [u]) which satisfy 

D y A - D X B + [A, M] + [L, B] = 0. (6) 

Whenever A and B are linearly independent, F is an immersion function for a 2D surface 
in the Lie algebra q. As proved in [H], three linearly independent terms which satisfy 
are 

d 

A = a—L + D X S + [S, L] + prv Q L G 0, (7) 

B = a-^M + D y S + [5, M] + prv Q M e 0, (8) 

where a = a(X) G C, S is an arbitrary g-valued function of A, y, x as well as the function 
u and its derivatives, and vq is a generalized symmetry of ([T|. Further, the q- valued 
function F can be explicitly integrated as 

F = a<$>- ld ^- + <$>- 1 S<$> + <$>- 1 prv Q <$>, (9) 

OA 

as long as vq is a generalized symmetry of the linear spectral problem, Q, as well as 
of the ODE §. 

In the next section, we consider a specific form for the ODE which includes 
differential equations for elliptic functions. We then present explicitly the Lax pair 
and find its wave function for the associated LSP. 

3. The Lax Pair and its wave function for a second-order ODE 

Consider a second-order differential equation given by 

u xx = \f'(u), f'(u) = A/(u) (10) 



for some function f'(u). It is straightforward to see that (10) admits the first integral 



{u x f = f(u), u x = t^J{u), e 2 = l, (11) 



and its solutions satisfy 
du 



x + x , x G E. (12) 



In the case that 



-}== = R(u, v^R) (13) 

vf( u ) 
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where R is a rational function of its arguments and P(u) is a polynomial of degree 3 or 



4, then the function u which solves (10) is the inverse of an elliptic integral [HI El El [16 
In particular, this is the case when f(u) is a polynomial of degree 3 or 4. 



In this section, we will construct a Lax pair for (10) in terms of matrix functions 
L, M taking values in sl(2, M) which satisfy the Lax equation ^ and find solutions for 
the wave functions $ G SX(2,R) which are solutions of the LSP Q. Let us make an 
assumption for the form of M 



M 



u x mi2 
u + X —u x 



det(M) = -g(X), 



(14) 



where g(X) is an arbitrary function of A. In what follows, we call g(X) the discriminate. 

(15) 



From (14), we obtain 

det(M) = -g(X) = -u\ - (u + X)m 12 , 

and making use of the first integral ( JlT| , we find that m\ 2 has the form 

f(u)-g(X) 



m 12 



u + X 



(16) 



Note that mi 2 , and hence M, are rational functions of A. Further, if f(u) is a polynomial 
in u, then m 12 will be a polynomial in u if and only if 



<?(A) = /(-A). 
Next, we solve for L written in the basis of sZ(2,R) 
L = l x e 2 + £ 2 e 2 + 4e 3 , 



(17) 



where £i i = 1, 2, 3 are some functions of A, x, u and u x . The matrices written in this 
basis are 



ei 














1 




1 


1 





, e 2 = 








, e 3 = 


-1 



Applying the total derivative D x to m 12 , we obtain 

/ /'(«) f(u)-g(X) 
D x {m l2 ) = — - + 

\ it + a (u + xy 

The Lax pair equation ^ becomes 
u^-iS-^-Uu + X) (-^ 



u. 



;is) 



(19) 



/(u)-g(A) 
«+A (m+A) 2 



M x + 2£ 2 M :C 



2£ 3 /(")~g( A ) 



u+X 



u x - 2t x u x + 2£ 3 (u + A) -u xx + ii 1 ^^ 1 + £ 2 (u + A) 







With the choice of functions, 

1 1 ff{u) 



2\u + X 



f(u)-g(X) 

(u + xy 



(20) 



(21) 
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(20) becomes 



u xx - \f{u) 










-u xx + \ f'{u) 











which is equivalent to the ODE (]10|). Thus, the matrices L, M of the form 

fjn)-g(X) 



u+X 



(«+A)2 

1 



M 



u + X 



u+X 

-u x 



(22) 



(23) 



satisfy the Lax equation (J2]). 

Next, we look for the most general solution of the wave function which satisfies the 
LSP Q. The components of the wave function $ are denoted by 



$ 



$n $i 2 
$21 $22 



G SL(2, 



with solutions of the LSP Q 

$11 = C101+ + C 2 0l_, 
$21 = Ci0 2 + + C 2 </>2-, 

and where 



$12 
$22 



C301+ + ^AVl-i 
C 3 4> 2 + + C 4 4>2- 



CiEK 

i = 1,2,3,4 



±\/gJX)+u x 



Vu + X 

2± = VM + A*±, 



± 



exp 



2/ 



2 (it + A) 



(24) 

(25) 
(26) 

(27) 
(28) 
(29) 



Here the choice of e comes from (jl 1 p . Note that, for the purposes of symmetry analysis 

th 
du 



it is sometimes useful to express the the integral appearing in ( 29 ) as 

dx 



edu 



2(u + X) J 2(u + X)u x J 2(« + A)v//g ^ 

by invoking (11). This integral will be explored further for particular choices of f(u) in 
the following sections. The derivatives of the functions \I/± have the simple form 



D y ^ ± = ±v^(A)^ ± , D x *± = ±; 



±- 



(31) 



~2{u + A) 

We will now check that the wave function $ as defined by (25 29) satisfies the LSP Q. 
In terms of the functions (p a ±, the differential equations Q to be verified are reduced 
to the following 

f'(u) f(u)-g(Xy 



D 
D 

Dy 

D 



xW2±i 



U + X 



(u + xy 



>2± 



0. 



''lit 



0. 



f(u)-g(X) ^ 
h± - u x (pi± H — 2 ± 



+ U 



X V2± 



u + X 
(u + A)0i ± 



0. 



0. 



(32) 
(33) 

(34) 
(35) 
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We verify (32) by computing the following total derivative 
DJ<t>i±) - 



u. 



1 (±y/g(X) + u x )u x ±y/g(X)(y/g(X) ± u x) 



yJu + X 2 yju + X(u + A) 



2(u + X)yju + X 



+ 



(u + X) 

' /'(") 

2{u + X)^ (u + xy 



— ^ Vu + A#± 
{u + X) 2 J 

-f(u)+g(X) 



>2±- 



Thus, we can directly observe that (32) is satisfied for the set of linearly independent 
solutions {0i±,02±}- Next, we verify (33) by computing the total derivative 



u x 



>2±, 



2\/u + X 

u x ±y / g(X) 
2Vu + X 



+ VuTx 



2(u + X) 



Thus, (|33|) holds. To verify (|34|), we note, using (|31|), that the dependence of $ on y is 



straightforward and so the computation is 



D 



U X 4>1± + 



U + X 



4>2± = (± Vflr(A) - u x )(f) 1± + 



<?(A) 



uz 



u + X 

<?(A) 



J 2± 



\Ju + X yu + X 



0. 



Thus, (34) holds. The last equation to verify is (|35j), 

Dy4> 2 ± + u x (j) 2 ± -{u+ A)0i± = (±a/#(A) + u x )(p 2 ± -{u+ A)0i± 

+t /~Tv\ _l ^ r~T\ (u + X)(±y^gJX~) +u x ) 

±W9W + u x)Vu + X 

\/u + X 



0. 



We have thus proven that the set of linearly independent solutions {<pi±,4>2±} satisfies 
(32 35). Finally, we note that the requirement det($) = 1 gives an algebraic constraint 
on the constants c,-, 



2^gJXj(c 2 c 3 - 



cic 4 ) 



In the previous section, we have assumed that the Lie algebra is sZ(2,R). In 
particular, we assumed that the dependent variable is a real function of a real variable, 
x, and the auxiliary variable, y, to be real as well. On the other hand, we can consider 
the complexification of the independent variable x, for example in the case of the V- 
Weierstrass function treated in section [6] which naturally extends to the complex plane. 
In this case, we can interpret the auxiliary variable y as the complex conjugate of x with 
the assumption that u(x) is analytic in x. The ODE (JT]) along with the requirement 
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that u be analytic gives a set of real PDE's for the real and imaginary parts of u. The 
above construction of L and M still hold, only now they take values in s/(2,C). The 
function / is assumed to be a complex function of u. Here the spectral parameter, A, 
takes values in C. 

For the purposes of constructing surfaces generated by the wave function $, we 
assume that it is a real function and x, y and A are also real. We choose the constants 
in $ to be 

1 1 

2' 



ci = c 2 



c 3 = -c 4 



so that the wave function $ simplifies to 



7=rr' 



01-) 
<t>2-) 



(36) 



and $ G SL(2, R) for all values of A G K. In particular, if g(X) < 0, then 



1,2 



so that their sum will be real and their difference, divided by ^g(X), will also be real. 



4. The induced surfaces 



For analytical descriptions of a 2D surfaces, the matrices A and B, given in equations 
([7]) and ^ are assumed to be linearly independent. With the wave function $, given 



by p4[ ) or (36), we can construct a surface F G s/(2,M). We construct separately three 
cases since the pair of matrices A and B corresponds to three types of symmetries. These 
are a conformal symmetry in the spectral parameter A (called the Sym-Tafel formula), 
gauge transformations of the wave function and the generalized symmetries of the ODE 



(10). These three types of symmetries yield different types surfaces. Let us consider 



each type of them individually. 



4-1- Sym- Tafel formula for immersion 

The first term in ^ and ^ corresponds to the Sym-Tafel formula for immersion which 
is given bv [2T| 122] 



(37) 



where a(A) is an arbitrary function of A. This surface F ST has tangent vectors of the 
form 



D X F 



D y F 



st $-i(_^_ M )$. 

OA 



(38) 



With L and M given by (23), the tangent vectors in (38) are linearly independent and 
so the function F ST gives an immersion of a 2D surface in the Lie algebra s/(2,R). 
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4-2. Surfaces associated with gauge symmetry 

The second term in ^ and ^ corresponds to the gauge symmetry of the linear spectral 
problem The surface F s corresponding to this gauge term can be integrated 

explicitly as QH [13] 

F s = $- 1 S(\,y,[u])$ (39) 

with tangent vectors 

D X F S = (D X S + [S, L}) $, D y F s = (D y S + [S, M]) $. (40) 

Again, for F to be an immersion, we require the linear independence of the tangent 
vectors. 

Note that, for any surface P e sZ(2,R), P can be expressed as 

p = Q^S® = F s (41) 
and hence F s represents a completely arbitrary surface immersed in the Lie algebra. 



For any sl(2, IR)-valued function 5*, written in the basis (18) as 

S = si(A,y, [u))ei + s 2 (X,y, [u])e 2 + s 3 (A,y, [u])e 3 , (42) 
the surface F s takes the form 



F s = si 



'ei - e 2 ) H 7 — , rr e 3 



4(« + A) V1 zy (« + A)V# 

+ S 2 { ^-^- + *-)% + " ^-) 2 ^ " - 

4 V 5 \fg J 



+ S 3 



+ (K. + v^)^+ + (u* - ej 



(tt _ <jr_) (( UiB + ^)^ + - ( Mz - v ^)$_) K + y/g)V\ - {u x + ^/g)^ 2 _ 

H -e 2 es 

2<? 2^ 



We can interpret the surface given in the form (41) as an arbitrary surface immersed in 
this Lie algebra written in the frame defined by conjugation by the wave function <E>, an 
element of the Lie group SX(2,R). 

4-3. Surfaces associated with generalized symmetries 

The third term of (J7J) and ^ is associated with generalized symmetries of equation 



(84). That is, suppose that there exists a generalized vector field with evolutionary 



representative 



v Q = Q[u}-^ (43) 



which is a generalized symmetry of (10) in the sense that 

prv Q [u xx -^f'(u)j =0, whenever u xx - ^ f'(u) = (44) 
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holds, where we have used the standard definition of the prolongation of a generalized 
vector field given as in the book by P. J. Olver [T9] , 



w t Q = Q ^ +DxQ _^ + DlQ JL + .... (45) 
The determining equations for Q are thus 

DlQ - \f"{u)Q = 0, whenever u xx - l -f{u) = 0. (46) 



For such a generalized symmetry (43), there exists a surface immersed in the Lie algebra 
q, say F®, with tangent vectors [H] 

D X F Q = $- 1 prtf Q L$, D y F Q = $- 1 prv Q M<$>. (47) 

Further, if the generalized symmetry is also a symmetry of the LSP in the sense that 

prv Q (D x <f> - L$) = 0, whenever D x f - L$ = (48) 
prv Q (D y § - Af$) = 0, whenever D y $ - M$ = 0, (49) 



then (47) can be integrated and the surface F® is given, up to a constant matrix of 



integration in sZ(2,M) by 

F Q = $-ipri? Q §. (50) 



The following characteristics Qj's are solutions of the determining equation (46) 

(i) Qx = u x 

(ii) Q 2 = u x J f(u)~zdu 

(iii) Q 3 = xu x + ju, this is only in the special case f(u) — c\ + c 2 u e for £ = 2(1 + I/7), 
7,ci,c 2 e E. 

We begin our consideration with Qi = u x . It can be observed directly that the 
differential equation (10) is invariant under translation in x. Thus, the vector field 

v Ux = u x — (51) 



is a symmetry of (10). In fact, for any G(y : A, [u]) which does not explicitly depend on 
x, the prolongation of v Ux acts as a total derivative, 

prv Um (G(y, A, [u]) = D x (G(y, A, [«])) ^ A, [«]) = 0. (52) 

Hence, we can see that 

prv Ux (u xx - ~/'(it)) = Dj; - = 0, whenever u xx - ^f'(u) = 0. (53) 

Similarly, the wave functions and the potential matrices L and M do not depend 
explicitly on x and so 

prv Ux (D x <& - L$) = D x (D x $ - L$) = 0, whenever D x $ - L$ = 0, (54) 

prv Ux (D y $ - M $) = A, (D^^ - M$) = 0, whenever - M$ = 0. (55) 
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Thus, for an arbitrary function f(u), the vector field v Ux is a symmetry of both (10) and 



the LSP Q and so according to [14], there exists a surface defined by the immersion 
function 



jpu x _ 

with tangent vectors 

D x F Ux = Q^D^L)® 
For the second case, we can verify that Q2 



D y F Ul 



(56) 

= &~ 1 D X (M)&. (57) 
w x / f(u)~^du solves the determining 



equation (|46j). The action of prvQ 2 on the LSP is given by 

-(*+ + ^_ 



prvQ 2 (D y $ - M$) 
prvQ 2 (D x $ - L&) 



It , 



0. 

■(*+ + tf- 





0(A)-i(¥+- 



^(A)-2(^ + -^_ 



(58) 



(59) 



2( M + A)iv^K 

Since these quantities do not vanish for all solutions of the LSP, the vector field vq 2 
is not a generalized symmetry of the LSP. Thus, while there exists an si (2, M)- valued 
immersion function F® 2 with tangent vectors 



DxF Q 2 = $-i prt ; Q2 (L)$, D y F Q2 = $->t;Q 2 (M)$, 



(60) 



the immersion function F® 2 is not of the form given in (50). 
Similarly, for the special case Q% 



xu x + 7« when f(u) = c\ + c^u , it is 



straightforward to verify that Q3 satisfies the determining equation (46). However, 
its action on the LSP gives 



prvQ 3 (D y <f> - M$) 
prvQ 3 (D x $ - L$) 



Ci(l + 7) 
\Ju + A 

Ci(l + 7) 



-(*+ + 




0(A)-*(¥- 



2(« + A)§ 







(61) 
(62) 



(*+ + 

Thus, unless ci = or 7 = — 1, vq 3 is not a symmetry of the LSP. In the former case, 
the differential equation reduces to the degenerate case u 2 x = c^ifi and in the latter case 
u is linear in x since u\ = c\ + c 2 . 

Thus, the surface associated with generalized vector field Vq 1 can be integrated 



explicitly and is given by (50) whereas the surfaces associated with the vq 2 and vq 3 are 



only of the form (50) for special cases. However, for all generalized symmetries there 



exists a surface with tangent vectors given by (47) and we can study their geometric 



properties based on the tangent vectors to the surface. To this end, in the next section 
we give a scalar product on the tangent spaces for the surfaces. 



4. 4- Induced metrics on the surfaces 



Here, we would like to introduce two possible choices for an induced metric on the 
tangents to the surface F G sZ(2,R). We take the basis for sZ(2,R) given by (18). A 
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first choice for a metric would be to decompose the matrix in the basis {ei, e^} and 
then to use the standard Euclidean metric. That is, given X, Y G sZ(2,IR) with 



X = X l a 



Y = Y l a 



1,2,3 



then the inner product and norm in Euclidean space is defined by 



(X, Y) = X'Y* 



\x\\ = y/xyx*. 



(63) 



(64) 



This constitutes an inner product on the tangent vectors for the surface F G sl(2, R). 

On the other hand, we can also use a symmetric bilinear product defined using 
the Killing form, B(X,Y). The main advantage of this form is that it is invariant 
under conjugation by the group and so, because of the form of the tangent vectors ([5]), 
the geometric quantities associated with the surfaces will be independent of the wave 
function $ G SL(2, R). The Killing form on sZ(2,R) is given, up to a normalization 
factor, by (see e.g. [T5] ) 

B(X,Y) = tr(XY). (65) 

In terms of the basis {ei,e2,es}, the matrices X, Y G si(2,R) and the Killing form 
B(X, Y) can be represented as the following 





X 1 




Y 1 




X = 


X 2 


, Y = 


Y 2 


B(X,Y) = X T BY, (66) 




X 3 




Y 3 





with 



B{X,Y) 



1 

1 
2 



(67) 



The Killing form has signature (2, 1) and so induces a pseudo-Euclidean metric on the 
tangents to the surface given by the immersion function, F G sl(2, R). The surface 
defined via the immersion function F with the Euclidean metric ( , ) is a Riemannian 
manifold while the surface with the Killing form B is a pseudo-Riemannian manifold 
0123]. 

As an example, for a general f(u) the first fundamental form for the surfaces F ST 
and F Ux with the pseudo-Euclidean metric are given by 



Ib(F st ) 
ff" 



f-g S'-sf 



and 

M \u + X {u + X) 2 

where for convenience, we denote 

/ = f(u), 



(u + X) 3 (u + X) 2 J 

Vf , V(f-g) 



I dxdy 



+2 



g' 



u + X) 3 

g = gW, f 



dxdy- 



r 

2 



v + X 



+ 



f-g 

(u + X) 2 

V(f-g) 



u + x (u + xy 



dy 2 , (68) 



dy 2 , (69) 



d -f(u), g' = ^g(X). 



du dX* 

In the Euclidean metric, the components of the metrics are significantly more 
complicated and can be found in Appendix A| Similarly, the components of the metric 
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Table 1. Constants in (70) for the Jacobian Elliptic functions. 



ki &2 Solutions of (70) 



1 

k' 2 
-k' 2 



-k 2 
k 2 
1 



sn(x, k) 
cn(x, k) 
dn(x, k) 



for F , in both the Euclidean and Pseudo-Euclidean metric, are too involved to write 



out in an illustrative fashion. However, they are directly computable from (43). In 
practice, the gauge term can be used to simplify the expressions for the surfaces. 



5. Jacobian elliptic functions and associated surfaces 

Consider the differential equation for the Jacobian Elliptic functions 



or alternatively 



-2k 2 u 3 + (k 2 - kx)u. 



(70) 



(71) 



Solutions of (|70|) are given by Jacobian elliptic functions, for the choice of constants 

(72) 



given in table [T] with 

k' 2 + k = l, 0<k,k'<l, 
In terms of the notation of section |3j the function f(u) takes the form 

f(u) = (l-u 2 )(k 1 + k 2 u 2 ). 

Choosing 

( 7 (A) = /(-A) = (l-A 2 )(A; 1 + M 2 ), 
the matrices L and M become 



M 



u x (u - X)(k 2 {u 2 + A 2 ) +ki — k 2 ) 
u + A — u x 

-3k 2 u 2 + 2Xk 2 u + h-k 2 - k 2 X 2 

1 



(73) 
(74) 

(75) 
(76) 



Taking the normalization of $ as in (|36|), the wave function takes the form 
$ = 



(y/gW-u x )^ + -(y/g(X)+u x )^- (y/g(X)+u x )^^-(y/g(X)-u x )9 + 



/u+A(* + +*_) 



2y/g(\)Vu+\ 
v/u+A(*_-* + ) 



where 



ty± = exp 



dx 



2{u + A) 



(77) 



(7f 
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Recall that, for the purpose of computing the action of the generalized vector field vq 



on $, it is convenient to use the identity (30) 

dx 



edu 



2{u + X) J 2 (u + \)y/(l-u^(k 1 + k 2 u 2 ) 



(79) 



when u is a solution of (|70|). The integrated form of (|79|) is 
/ 2(u 



edu 



>(u + \)y/(l - u 2 )^ + k 2 u 2 ) 2\^k 
-.tank 




. _! / (k 2 - k 2 - 2k 2 X 2 )u 2 + (k 2 - h)\ 2 + 2h 



A^gJX) \ 2^gJX}^/(l-u^(k 1 + k 2 u 2 ) 

where II (u, a, b) is the normal elliptic integral of the third kind, see e.g. [6] 

II (u, a 2 , k) = [ r 

V 1 J (1 -a 2 t 2 )vT^tVl - kH 2 

The functions ty± are then given by 



+ c , (80) 



(81) 



ty± = exp 



±v/9(A)|!,+ vsr n t"'^v^ ?l+c ° 



X 



2^gJX)^{l-u 2 ){k 1 + k 2 u 2 ) + (k 2 -k 2 - 2k 2 X 2 )u 2 + (k 2 - k^X 2 + 2kt 



82) 



_2 v ^7(A)v / ( 1 -u 2 )(h + k 2 u 2 ) - (k 2 -k 2 - 2k 2 \ 2 )u 2 - (k 2 - h)X 2 - 2h 

In the graphs below, we choose the integration constant Cq so that ^±(0, 0) = 1. 

From the wave function $, it is immediate to compute the analytical form of the 
surface generated by the terms described in section |4j 

F = (a(X)^ 1 ^ + Q^SQ + bQ^D^j e sl(2, E) a(A), b e E. 

Below, we have presented several graphs of the induced surfaces plotted with the help 
of MAPLE. The surfaces in figure 1 have constants chosen so that the discriminant is 
negative, g(X) < 0, and so the surfaces behave like trigonometric functions whereas the 
surfaces in figure 2 have positive discriminant, g(X) > 0, and so have exponential type 
behavior. 

In the metric induced by the Killing form, the first fundamental form for the surface 
F ST is 



Ib{F st ) 
and fc 
Ib(F u * 



^(k 2 (u - X)dxdy - [k 2 (u 2 



2Xu + 3A 2 - 2) + l] dy 2 . 



and for F Ux is 

k 2 b 2 



+ - 



(3m - A)(l - u 2 )(k 2 u 2 + k' 2 )dxdy 
k 4 u 6 + 2u 5 k 4 X - u 4 k 4 X 2 - {Ak A X - 2Xk 2 ) u 3 + (3A; 2 + 2A; 4 A 2 - 3k 4 - k 2 X 2 ) 



u 



- (2Xk 2 - 2k 4 X) u - k 4 X 2 + 2k 4 + k 2 X 2 -3k 2 + l 



dy 2 
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c3 




e3 




F ST : A = 1.2, k = 0, F ST : A = 1.2, fc = 0.5, F ST : A = 1.2, k = 0.8, 




e3 01 




e3 




F Ux : A = 1.2, fc = 0, F Ux : A = 1.2, jfe = 0.5, : A = 1.2, k = 0.8, 



Figure 1. Surfaces i 7,52 " and for u — sn(x, k) with the discriminant g(X) < and 
x and y € [—8, 8]. The axes indicate the components of the immersion function in the 



basis ( 18 1 




F ST : A = 0.5, fc = 0.2 



F Ux : A = 0.5, fc = 0.2 



Figure 2. Surfaces F ST and i 71 "^ for u — sn(x,k) with the discriminant g(X) > 0, 
x E [—20,20] and y E [—5,5]. The axes indicate the components of the immersion 



function in the basis ( 18 ) 
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We can also study the surfaces associated with generalized symmetries Q 2 = 
u x J f(u)~zdu and Q 3 = xu x + 7W, in the special case when k\ = —k 2 , 7 = 1 and 



f(u) = k\(l — u 4 ), but the immersion functions are not of the form (50) since the vector 
fields vq 2 and vq 3 are not generalized symmetries of the LSP. However, we can study 
the geometric properties of these surfaces via their tangent vectors. 



6. P-Weierstrass elliptic function and associated surfaces 

Consider the differential equation for the P-Weierstrass elliptic function 

u x = -\/4m 3 - g 2 u - 03, 
or alternatively 

u xx - 6w 2 + y = 0. 
In terms of the notation of section |3j 

f(u) = 4m 3 - g 2 u - g 3 . 

If we choose 

0(A) = /(-A), 
the Lax pair L, M becomes 

// , -4m 2 + 4Aw - 4A 2 + g 2 



M 



1 



u + A 

4m-2A 

1 



-u. 



(83) 
(84) 
(85) 
(86) 

(87) 



Taking the normalization of $ as in (36), the wave function takes the form given 



by (|77|) but differs from the wave functions from the previous section since the functions 

dx 



^± are instead given by 



ty± = exp 



9) 



2(u + A) 

Again, for computing the action of the generalized vector vq field on $, it is convenient 
to use the identity 

dx f du 



2{u + A) 



2(u + X) v /Au i -g 2 u + g 3 ' 



(90) 



when u is a solution of (84). The integrated form of (90) is 
du 



2(u + A) a/ 4m 3 - g 2 u - g 3 
-1 



2(A-ai)(2ai + a 2 ) 



n 



u + ai a± — a 2 
ai — a 2 A — a\ ' V 2ai + a 2 



ai — a 2 



+ c , 



(91) 



Surfaces immersed in Lie algebras associated with elliptic integrals 



17 



where, for simplicity, we have used the constants ai, a 2 to represent the roots of the 
polynomial f(u) as in 

f(u) = 4w 3 - g 2 u - g 3 = 4(w + a x ){u + a 2 ){u -a x - a 2 ). (92) 

Here, II is the solution of the elliptic integral defined by ( fSlp . In the graphs below, we 
choose the integration constant Cq so that ^±(0.5,0) = 1. 

From the wave function, <3> it is immediate to compute the analytical form of the 
surface generated by the terms described in section |4j 

F = ^a(A)$" 1 ^$ + Q^S® + b^D^j e sl(2, E) a(A), feel. 

In figure 3, we have included graphs of the induced surfaces for P-Weierstrass functions 
with g 2 — 0, #3 = 1, i.e. on an equianharmonic lattice [3]. Again, the surfaces exhibit 
periodic behavior when g(X) < and exponential behavior when g(X) > 0. We can 
also study the surfaces associated with generalized symmetries Q 2 = u x f f{u)~^du and 
Qs = xu x + 7«, in the special case where g 2 = 0, 7 = 2 and f{u) = 4m 3 — g 3 . However, 



the immersion functions are not of the form (50) except in the latter case when g 3 is 



also equal to zero and the solutions of (84) are rational functions of x. 

In the metric induced by the Killing form, the first fundamental form for the surface 
F ST is 

I B (F ST ) = -2a 2 (A) (dxdy + 3a 2 (2A - u)dy 2 ) . 
and, for F Ux , is 



I B {F Ux ) = 46 2 



(4m 3 - g 2 u - g 3 )dxdy 



+ (^2n 4 + 8Xu 3 + g 2 u 2 - 2{g 2 \ - 2g 3 )u - 2Xg 3 + -g 2 2 j dy 

Both fundamental forms depend on the solution u and the spectral parameter A. Note 
that, in both cases, the tangent vectors in the x-direction are null vectors in the pseudo- 
Euclidean metric. 



7. Final remarks 



In this paper, we have discussed certain classes of surfaces immersed in Lie algebras 
associated with elliptic integrals. This problem has been studied recently for PDE's 
by the authors in [H] using a symmetry characterization of continuous deformations of 
soliton surfaces immersed in Lie algebras based on the formalism of generalized vector 
fields and their prolongation structure. The necessary and sufficient conditions for the 
existence of such surfaces in terms of the invariance conditions has been established 
for integrable PDE's. In this context, we have adapted the proposed procedure for 
integrable ODE's admitting a Lax representation ^ and shown, as in the PDE case, 
that the problem requires the examination of conformal symmetries in the spectral 
parameter, gauge transformations of the wave function and generalized symmetries 
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e3 




F ST : A = 1, g(X) 





F Ux ■ A = 1, g{\) 



F Ux : A = -5, g(X) = 31 



Figure 3. Surfaces F u * for u = T(x,0,l) with x E [0.2,3] and y E [-71-/5(1), n/g(l)} 
for the negative discriminate cases and x E [—1,1] V E [—0.5,0.5] for the positive 
discriminant cases. The axes indicate the components of the immersion function in 



the basis ( 18 ) 



of the associated model and its LSP. To perform this symmetry analysis, we have 
constructed a Lax pair for a second-order ODE which includes, among others, the case 
of elliptic integrals. Next, we solved the LSP and found explicitly the most general form 
of the wave function. Next, we constructed surfaces in s/(2,R) by analytic methods 



for a general case where the right hand side of (11) is an arbitrary function f(u). We 
were able to explicitly integrate these ODE's in terms of elliptic integrals if f(u) is a 
polynomial in u of degree 3 or 4 and give explicit forms of the corresponding soliton 
surfaces for the Jacobian and P-Weierstrass elliptic functions. A geometrical analysis of 
these surfaces has been performed by using an appropriate inner product which allows 
for the construction of Riemannian and pseudo-Riemannian manifolds. The elaborate 
procedure was applied to examples and we have given the first fundamental form for the 
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surfaces as well as graphs of the surfaces for a range of parameters leading to diverse 
types of surfaces. 

Additional questions which could be asked involve the geometric properties of the 
surfaces associated with elliptic integrals. It would be especially interesting to determine 
some global characteristics of the soliton surfaces defined by the immersion function F. 
These quantities, such as the Willmore functional or Euler-Poincare characteristics as 
well as the Gaussian and mean curvatures, can be calculated but the expressions are 
rather involved, so we omit them here. It is an open question as to whether, by judicious 
choice of the linearly independent terms in the matrices A and B given by ^ and (|8|, 
the equations can be simplified. Finally, another interesting avenue for future research 
could include the application of the methods presented here to other integrable equations 
ODE's, for example hyperelliptic functions, which admit a Lax representation and to 
study their soliton surfaces. 
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Appendix A. Components of the first fundamental form for surfaces in the 
Euclidean metric 

The components of the first fundamental form for the surface induced via the Sym-Tafel 
formula for immersion, in the Euclidean metric, are given by 



(f- g ')(u + X)-2(f-g) 



64(u + X) 2 g(X) 

x ((g 2 + g + + **_) + 4(g 2 - + vfr 2 J + 2(3g 2 -g + 3)), 

2(3g 2 -g + 3)(u + \)g' - Ag(g 2 + g + 1) 



,st zpST\ 2(f-g)-(f'-g')(u + X) 
" ' 32(u + Xfg{X) 



-2 (* 4 + - 1) ((g 2 +g + 1)(1 + *i) + 2{g 2 - 1)^ 2 _) u x 

+(g 2 + g + l)(g'(u + A) - 2g){^% + + % 2 - l)(g'(u + A) - g)^ 2 + + * 2 _) 
<^f T > F v T ) = ~ A ( 5 ~ l L ((9 2 + 9 + Win + A) - 2g)(l + *i) + 2g'{g 2 -l)(u + \)* 2 _) 

(g 2 + g + l)((g'(u + A) - 2g) 2 + Agf) (^ + + * 4 _) 



u, 



16(u + X) 2 g 2 



X 



+4g'(g 2 - l)(u + X)(g'(u + A) - 2g)M + *i) 



(A.l) 
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+2(u + A) 2 (3 -g + 3g 2 )g /2 - 8g(g 2 -g + l)(u + \)g' 
where for convenience, we denote 

'"I'M- *' = ^ A »- 

Similarly, the components of the first fundamental form for the surface induced via 
generalized symmetries are 

(f q pQ) = Q 2 (f"(u + X) 2 -2(u + X)f + 2(f-g)) 2 
[ x ' x ' ~ Q4(u + \) A g 2 

x [(g 2 + g + 1)(* 4 + + * 4 _) + Af(g 2 - 1)(* 2 + + %2 _ ff + 3) ] > 
(f q F Q\ = Q{f\u + \) 2 -2{u + \)f + 2{f-g)) 2 ) 
K x ' » ; " 32(u + A)V 

x Qg(g 2 + 1) + 2y/g(D x Q - Q Ma; ) ((/ + g + l)(\I>t " + - !)(*+ - 
+ (2D x Qu x - Qf\u + A) + 2gQ) 

x [fa 2 + g + 1)(M/ 4 + H/ 4 „) + 2(<? 2 - 1)(^ + M/ 2 ) + 2(3/ - g + 3)] 

Q Q = (g 2 + g + l) 
K v ' v 1 \Qg\u + \Y 

x [4v^(/'(« + A) + 2g)u x Q 2 + 2(£> 2 QK - 2(f(u + X) + 2g + 2f)Q x Q] (^ 4 + - ^ 4 _) 

+ [((/'(« + A) + 2^) 2 + 4<//) Q 2 - 4(/'(« + A) + Ag)D x QQu x + 4(/ + sX^Q) 2 ] (M> 4 + + M> 4 _) 

- (4(/'(u + A) + g)D x QQu x + (f'(u + A) + 2< ? )/'Q 2 + (« + \)f{D x Q) 2 ) (V 2 + + ^ 2 ) 
(3g 2 -g + 3)(f'D x QQu x + f(D x Q) 2 ) + ^ - ff + 1XA.Q) 2 



S 2 



(3<? 2 -<? + 3)/' 2 , (,/-• // • i )/' , (/r // • Hi// f)\ 



+ V 8^ + 2,7(11 + A) + 2#(u + A) 2 J Q - (A ' 2) 

For the vector field v Ux the components simplify to 

(f q F Q\ = f(f"(u + \) 2 -2( U + \)f> + 2(f-g)) 2 
K x ' " ; " 64(« + A) V 

x [(g 2 + g + l)(«t + + W " + + % 2 " 5 + 3)] , 
(F u x F u x) = U"{u + \) 2 -2{u + \)r + 2{f-g)) 2 ) 
" 32f?§(n + A) 3 

x (/'(« + A) - 2f)u x {{g 2 + g + 1)(* 4 + - * 4 _) + 2(</ 2 - 1)(* 2 - * 2 )) 

2^/ ((/ + 5 + l)(n + * 4 -) + 2(/ - ml + * 2 ) + 2(g 2 -g + 1)) 
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1 



k a + g + D(2/ -/'(« + A) 

(/> + A) 2 - 4/'/(« + A) + 4/(/ + + <? + l)(*t + * 4 _) 



(« + A) 2 ^L 

-2(/> + A) 2 - Aff(u + A) + 4/(/ - g))(g 2 - g + 1) 



(A.3) 
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